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ABSTRACT 

We present observations of the A7IV-V star Alderamin (a Cep, HR 8162, 
HD 203280) from the Georgia State University CHARA Array. These infrared 
interferometric angular size measurements indicate a non-circular projected disk 
brightness distribution for this known rapid rotator. The interferometric obser- 
vations are modeled as arising from an elongated rigid atmosphere, with appar- 
ent polar and equatorial radii of r p = 0.67531q;q^5 milliarcseconds (mas) and 
r e = 0.8767^o'oi83 mas, respectively, for a difference of 201 ± 32 microarcseconds 
(/zas), and with an axial ratio of r e /r p = 1.298 ± 0.051. Using the Hipparcos 
distance of 14.96 ±0.11 pc, these angular measures translate to 2.18 ± 0.05 and 
2.82 ± 0.10 R Q . The inclination of Alderamin to the line of sight indicated by 
this modeling is effectively edge-on (i = 88.2+J 3 8 3 ). The star has a true rotational 
velocity of 283 ± 10 km/s (~ 83% of breakup velocity), and a polar temperature 
of roughly 8400 K. Significantly, a necessary aspect of this modeling is a determi- 
nation of the gravity darkening coefficient, which at a value of (3 = 0.084+q 049 is 
consistent with a convective photosphere, as expected for an A7IV-V star. Our 
detailed characterization of this object allows us to investigate various scenar- 
ios for the angular momentum history of Alderamin and the appropriateness of 
certain stellar evolution models. 

Subject headings: stars: individual: Alderamin, infrared: stars, stars: fundamen- 
tal parameters, techniques: interferometric 

1. Introduction 

The Georgia State University's (GSU) Center for High Angular Resolution Astronomy 
(CHARA) Array is a six-element optical/infrared interferometer located on Mt. Wilson in 



1 For reprints, please contact: gerard@ipac.caltech.edu. 
2 Michelson Fellow. 
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southern California. The CHARA Array has six 1-m telescopes operational and recently 
completed its first full year of science observations, including observations of stellar diame- 
ters, young stellar objects, and rapidly rotating stars. A companion paper (ten Brummelaar 
et al. 2005) describes the full compliment of technical details of the instrument, and McAl- 
ister et al. (2005) details the first science results from the instrument, on the rapid rotator 
Regulus. CHARA's operational status of having the longest H- & if-band baselines in the 
world make it uniquely well suited for observations of main sequence star absolute diameters. 

The star Alderamin (a Cep, HR 8162, HD 203280) is a well-studied object, being one 
of the 20 brightest A-type stars in the sky (Hofneit & Warren 1995) and one of the nearer 
stars to the Sun (Perryman et al. 1997; Cox 2000). Originally classified as an A2n (Douglas 
1926), it is now accepted to be an A7IV-V main sequence star (Johnson & Morgan 1953) and 
has been known to be a rapid rotator for over 50 years (Slettebak 1955). The measurements 
of the star's apparent rotational velocity (v sin?) range from 180-200 km s^ 1 (Gray 1980; 
Abt & Morrell 1995) up to 245-265 km s" 1 (Bernacca & Perinotto 1970; Abt & Moyd 
1973), depending upon the spectral lines used in the investigation. These values of wsini 
are a substantial fraction of an A7V star's critical velocity of no more than 415 km s _1 (a 
velocity which decreases as the object evolves and increases in size) (Slettebak 1966), where 
centripetal acceleration at the stellar equator equals gravitational acceleration. 

Recent advances in interferometric observations have allowed for direct observation of 
rotationally oblate main sequence stellar surfaces, including the first such observation by 
van Belle et al. (2001) of the A7IV-V rapid rotator Altair, confirmation of that object's 
latitude-dependent temperature structure (Ohishi et al. 2004; Peterson et al. 2005), and the 
CHARA Array observations of Regulus (McAlister et al. 2005). Similar observations of the 
related Be star phenomenon have been also achieved by interferometers (Quirrenbach et al. 
1994; Domiciano de Souza et al. 2003). 

The effects of stellar rotation have been measured spectroscopically for almost a century, 
beginning with Schlesinger (1909, 1911). Models of rotating stars have explored the impact of 
rotation upon both stellar effective temperature (Slettebak 1949) and stellar shape (Collins 
1963, 1965; Collins & Harrington 1966). Recently, models have begun to incorporate the 
effects of differential rotation as a function of stellar latitude (Zahn 1992). Rotation impacts 
important observable parameters such as photometry (Collins & Smith 1985) and surface 
brightness distributions as originally shown by von Zeipel (1924a,b). Rotation has non-trivial 
implications upon stellar evolution, as explored in the various papers by, among others, Claret 
and Maeder (cf. Martin & Claret (1996); Claret (2000); Maeder (1997, 2000)). 

Herein we report the determination of the overall diameter and projected shape of 
Alderamin upon the sky from near-infrared, long-baseline interferometric measurements 
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taken with the CHARA Array. Direct observation of the stellar disk can provide unique 
insight into basic stellar parameters. The measured angular size in conjunction with the 
bolometric flux and distance yields constraints on parameters such as latitude-dependent 
local effective temperature 1 and linear radii, both of which remain quantities poorly estab- 
lished empirically for virtually all stars. Upon fitting a family of rotating models for the 
projected stellar photosphere upon the sky, we further demonstrate that a unique value for 
v sin % may be derived from the interferometric data. 

The CHARA Array observations that produced these results are discussed in §2, detail- 
ing source selection and observation. In §3, we detail supporting spectral energy distribution 
fits which constrain stellar parameters appropriate for this analysis. The circular symmetry 
of our check star is then established in §4, which also allow us to characterize the uncertainties 
inherent in the visibility data. In §5, Alderamin's departure from a circular on-sky bright- 
ness distribution is established, and in §6, we demonstrate that inclination, true rotational 
velocity, and other astrophysical parameters may be derived from Alderamin's oblateness by 
fitting the data with the appropriate family of Roche models. Finally, in §7, we examine 
the astrophysical implications of the best-fit model, such as the possible angular momentum 
history of the object. 



The interferometric observable used for these measurements is the fringe contrast or 
visibility (squared) of an observed brightness distribution on the sky. Normalized in the 
interval [0 : 1], a uniform disk single star exhibits monochromatic visibility modulus in a 
uniform disk model given by 



where J\ is the first-order Bessel function, B is the projected baseline vector magnitude at 
the star position, 9ud is the apparent uniform disk angular diameter of the star, and A is 
the wavelength of the interferometric observation. 

Alderamin was observed in the K s -b&nd (0.30/im wide centered at 2.15/xm) contem- 
poraneously with a primary calibration star, HD197373, and a check star, HD 211833, by 



1 We note that a star's effective temperature is a globally denned quantity in terms of the stellar luminosity, 
L = 4tt<tR 2 T fff . For denoting the temperature associated with specific surface elements of the stellar 
surface, we will follow the convention of Collins (1963, 1965) and use the term local effective temperature. 
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the CHARA Array on 8 nights between 2004 Jun 17 and 2004 Jun 29. Observations of 
Alderamin were always bracketed within ~ 20 minutes with the calibration source, and ev- 
ery other Alderamin-calibration set included an observation of HD2 11833. On the nights 
of Jun 17, 19, 20, and 21, Alderamin was observed with the Wl and El pair of CHARA 
telescopes; on Jun 26-29, the SI and El pair was utilized; details are given in Table 1. Data 
collection on the nights between Jun 22 and 25 was attempted with the Wl and SI pairing, 
but due to weather and instrumental difficulties, these data were of insufficient quality for 
this study and were discarded. Other nights between Jun 17 and Jun 29 were lost to weather 
and instrument problems. Our check star, HD 211833, was selected on the basis of an ex- 
pected angular size similar to Alderamin, in addition to an expectation of circular symmetry 
in its sky projection, based upon its low vsini. Alderamin and HD 211833, along with an 
"unresolved" calibration object, HD 197373, were observed multiple times during each of 
these nights, and each observation set, or scan, was approximately 180 s long, consisting of 
200 scans. Following commonly accepted & expected optical interferometry practice (van 
Belle & van Belle 2005), we attempted to achieve absolute instrument calibration through use 
of a calibration object that matched the instrument's resolution limit and limiting accuracy, 
as demonstrated by night-to-night repeatability. 

For a 0.41 mas calibration source such as HD 197373, the CHARA Array's longest base- 
line at 330m should give a raw V 2 of 0.80 before instrumental and atmospheric degradation. 
For each scan we computed a mean V 2 -value from the scan data, and the error in the V 2 
estimate from the rms internal scatter (ten Brummelaar et al. 2005). Alderamin was always 
observed in combination with its calibration or check sources within 5 deg and 6.8 deg on the 
sky, respectfully. The calibration source HD 197373 is expected to be nearly unresolved by 
the interferometer with a predicted angular size of 0.41±0.04 mas (computed in §3); expected 
angular size and error were based upon fitting template spectral energy distributions of the 
proper spectral type from Pickles (1998) to available broadband photometry, particularly in 
the near-infrared (Gezari et al. 1996; Cutri et al. 2003). Since many stars deviate signifi- 

Table 1. CHARA baselines utilized for observing Alderamin and associated sources. 



Baseline 


Projected baselines 


Projected sky angles a 


Dates 


Wl-El 


279-312m 


70-120° 


15-21 Jun 2004 


Sl-El 


250-304m 


0-75° 


26-29 Jun 2004 



PA is east of north. 
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cantly from blackbody behavior and/or have significant reddening (Blackwell & Lynas-Gray 
1998), we expect this approach to provide significantly better estimates of calibrator angular 
size than a simple blackbody fit. These objects were additionally selected to be slow appar- 
ent rotators, with vsini < 30 km s _1 (Uesugi & Fukuda 1982; Henry et al. 2000). Table 2 
summarizes the general parameters for the objects observed in this investigation. 

Interferometer U 2 's were obtained by recording a photometric signal of the two telescope 
combined beam as the interferometer delay lines were slewed through the white light fringe 
position on the sky at a pre-set group velocity. This signal was normalized with a low 
pass filter, and a power spectrum was computed. Using calibration scans of the individual 
beams and the closed shutter, the noise bias was removed from the power spectrum, and the 
integration of the power results in an estimate of instrumental V 2 . An additional correction 
for atmospheric turbulence visibility bias is also used to adjust this V 2 ; these steps are 
all described in mathematical detail in ten Brummelaar et al. (2005). Formal errors for 
each measure of V 2 were established from the measurement statistics associated with these 
photometric signals. 

The atmospheric calibration of Alderamin V 2 data on the sky is performed by estimating 
the interferometer system visibility (V 2 ys ) using the calibration source with model angular 
diameters and then normalizing the raw Alderamin visibility by V 2 to estimate the V 2 mea- 
sured by an ideal interferometer at that epoch (Mozurkewich et al. 1991; Boden et al. 1998). 
Multiple observations of the calibration source were averaged together in a time-weighted 
sense, with the error variance being doubled for a one hour time separation. Uncertainties 
in the system visibility and the calibrated target visibility were propagated though the data 
stream using standard error-propagation calculations. This atmospheric calibration process 
was accomplished through use of the publicly available wbCalib program 2 . 

The formal errors reported by the instrument measurement process are always smaller 
than those created by the atmospheric calibration process, which will be discussed further in 
§4. More detail on the CHARA Array's target and calibrator selection, data reduction and 
technical aspects is available in the literature (McAlister et al. 2005; ten Brummelaar et al. 
2005). Calibrating our Alderamin data set with respect to the calibration object listed in 
Table 2 results in a total of 41 calibrated scans on Alderamin over the 8 observing nights in 
2004, and 22 calibrated scans on HD 211833 over the same nights. Our calibrated Alderamin 
V 2 measurements are summarized in Table 3. 



2 Dctailcd documentation and downloads available online at http://msc.caltech.edu. 
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3. Spectral Energy Distribution Fitting 

For each of the three stars observed in this investigation, a spectral energy distribution 
(SED) fit was performed. This fit was accomplished using photometry available in the liter- 
ature as the input values, with template spectra appropriate for the spectral types indicated 
for the stars in question. The template spectra, from Pickles (1998), were adjusted by the fit- 
ting routine to account for overall flux level, wavelength-dependent reddening, and expected 
angular size. Reddening corrections were based upon the empirical reddening determination 
described by Mathis (1980), which differs little from van de Hulst's theoretical reddening 
curve number 15 (Johnson 1968; Dyck et al. 1996). Both narrowband and wideband pho- 
tometry in the 0.3 /nn to 3 /zm were used as available, including Johnson UBV (Eggen 1963, 
1972; Moreno 1971; Oja 1996), Stromgren ubvy(3 (Crawford et al. 1966; Piirola 1976), 2Mass 
JHK S (Cutri et al. 2003), Vilnius UPXYZS (Zdanavicius et al. 1969, 1972), and WBVR 
(Kornilov et al. 1991); flux calibrations were based upon the values given in Cox (2000). The 
results of the fitting are given in Table 5, and an example SED fitting plot is given in Figure 
1. 

The utility of this fitting was twofold. First, for our calibration source, HD 197373, an a 
priori estimate of its size is necessary to account for residual resolution that may be afforded 
by the interferometer's extraordinarily long baselines. With an expected limb darkened size 
of 9 EST = 0.412 ±0.019 from the SED fit, HD 197373 has a predicted V 2 of 79.8 ± 1.7% for a 
330-m baseline used at 2.2 /im; we shall consider this size effectively identical to its uniform 
disk size, since for a F-type size, the difference between the two is at the 1% level (Davis et 
al. 2000; Claret & Hauschildt 2003), which is far less than our size estimate error. Ideally, 
a calibration source would be sufficiently point-like that its V 2 would be indistinguishable 
from unity, but unfortunately the current system sensitivity does not afford that option. 



Table 2. Stars observed with CHARA. 



Source 


Q a 
VSED 


Distance from 


Spectral 


vsini 


Notes 




(mas) 


Alderamin (deg) 


Type 


(km s _1 ) 




Alderamin 


1.36 ±0.04 




A7IV-V 


~200 


Primary target 


HD 197373 


0.412 ±0.019 


5.0 


F6IV 


30 


Primary calibrator 


HD 211833 


1.34 ±0.06 


6.8 


K3III 


2 


Resolved check star 



a Estimated angular size from SED fitting as described in §3. 
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Table 3. The observed data for Alderamin. 



Projected Position Hour Normalized Uniform Disk 
MJD Baseline Angle Angle V 2 Ang. Size 







(m) 


(deg) a 


(hr) 




b 






(mas) 




53173. 


.363 


279. 


.16 


112.8 


-2.73 





.0323 ± 0. 


.0068 


1 


.579 ± 


.033 


53173. 


.385 


287. 


.16 


105.3 


-2.21 


0. 


.0145 ± 0. 


.0034 


1 


.632 ± 


.027 


53173. 


.415 


296. 


88 


95.3 


-1.48 


0. 


0062 ± 


.0015 


1 


.651 ±0 


.018 


53173. 


448 


305. 


.03 


84.9 


-0.70 


0. 


.0120 ±0. 


.0023 


1 


.554 ± 


.018 


53173. 


472 


309. 


.28 


77.4 


-0.13 





.0045 ± 0. 


.0009 


1 


.605 ± 


,013 


53175 


376 


285. 


87 


106.5 


-2.29 


0. 


0167 ±0. 


.0038 


1 


.624 ± 


.027 


53175. 


.421 


299. 


.99 


91.7 


-1.21 


0. 


.0088 ± 0. 


.0020 


1 


.607 ± 


.020 


53175. 


.434 


303. 


,27 


87.5 


-0.89 


0. 


0073 ± 


.0017 


1 


.604 ± 


.019 


53175. 


.473 


310. 


.25 


75.1 


0.05 


0. 


0059 ± 


0013 


1 


.583 ± 


.016 


53176. 


383 


289. 


.27 


103.2 


-2.06 





.0124 ±0. 


.0034 


1 


.635 ± 


.029 


53176. 


.398 


294. 


.11 


98.3 


-1.70 


0. 


.0108 ±0. 


.0030 


1 


.622 ± 


.027 


53176 


423 


301. 


.20 


90.2 


-1.10 


0. 


0082 ±0 


0023 


1 


.606 ± 


.024 


53176. 


.437 


304. 


.40 


85.8 


-0.77 


0. 


.0051 ±0. 


.0014 


1 


.623 ±0 


.020 


53176. 


.462 


309. 


.13 


77.7 


-0.15 


0. 


.0045 ± 0. 


.0012 


1 


.606 ± 


.018 


53176 


.490 


312. 


22 


69.0 


0.51 


0. 


0029 ± 


.0008 


1 


.613 ±0 


.014 


53177. 


437 


305. 


.05 


84.9 


-0.70 


0. 


.0049 ± 0. 


.0067 


1 


.621 ±0 


.138 


53177. 


456 


308. 


.58 


78.8 


-0.24 


0. 


.0064 ± 0. 


.0057 


1 


.586 ± 


.099 


53182 


348 


269. 


01 


50.4 


-2.51 


0. 


0425 ± 


0168 


1 


.594 ± 


.077 


53182. 


446 


298. 


31 


26.4 


-0.16 


0. 


0354 ± 


0146 


1 


.464 ± 


.068 


53182. 


.466 


301. 


.61 


21.3 


0.33 


0. 


.0690 ± 0. 


.0269 


1 


.341 ±0 


.082 


53182 


482 


303. 


69 


17.2 


0.72 


0. 


,0495 ± 


.0195 


1 


.388 ± 


.072 


53183. 


.306 


249. 


91 


60.1 


-3.45 


0. 


1114 ±0 


.0247 


1 


.502 ± 


.064 


53183. 


.322 


258. 


.10 


56.2 


-3.07 





.0976 ± 0. 


.0216 


1 


.488 ± 


.058 


53183 


352 


271. 


72 


48.8 


-2.35 


o 


0799 4- n 


0150 


\ 


.479 i 


.046 


53183 


368 


277. 


82 


45.0 


-1.97 


0. 


0607 ±0. 


0128 


1 


.481 ± 


.042 


53184. 


327 


261. 


.74 


54.3 


-2.89 


0. 


.0993 ± 


.0173 


1 


.463 ± 


.045 


53184. 


341 


268. 


.24 


50.8 


-2.55 


0. 


.0806 ± 


.0141 


1 


.475 ± 


.040 


53184 


.358 


275. 


.42 


46.5 


-2.13 


0. 


0660 ± 


.0121 


1 


.478 ± 


.038 


53184. 


381 


283. 


.38 


41.1 


-1.59 





0761 ±0. 


.0129 


1 


.408 ± 


.036 


53184. 


397 


288. 


.19 


37.2 


-1.21 


0. 


.0707 ±0. 


.0127 


1 


.399 ± 


.037 


53184. 


421 


294. 


30 


31.3 


-0.63 


0. 


0608 ± 


.0136 


1 


.398 ± 


.043 


53184. 


433 


297. 


02 


28.1 


-0.32 


0. 


0664 ± 


0134 


1 


.369 ± 


.039 


53184. 


.458 


301. 


.24 


22.0 


0.27 


0. 


.0715 ±0. 


.0133 


1 


.337 ± 


.037 


53185. 


375 


282. 


.19 


41.9 


-1.67 


0. 


0874 ± 0. 


.0231 


1 


.385 ± 


.062 


53185. 


388 


286. 


.43 


38.7 


-1.35 


0. 


0467 ±0. 


.0142 


1 


.482 ± 


.056 


53185. 


.410 


292. 


.33 


33.3 


-0.83 





.0725 ± 0. 


.0211 


1 


.374 ± 


.062 


53185 


419 


294. 


64 


30.9 


-0.59 


0. 


0702 ± 


0213 


1 


.370 ± 


.063 


53185 


437 


298. 


14 


26.6 


-0.18 


0. 


0770 ± 


0235 


1 


.336 ± 


.065 


53185. 


.448 


300. 


.09 


23.8 


0.09 


0. 


.0550 ± 0. 


.0159 


1 


.388 ± 


.053 


53185. 


.469 


303. 


.05 


18.6 


0.59 


0. 


.0398 ± 


.0121 


1 


.425 ± 


.049 


53185. 


.481 


304. 


40 


15.5 


0.88 


0. 


0505 ± 0. 


.0151 


1 


.382 ± 


.053 



a PA is east of north. 

b Errors have been normalized as discussed in §4. 
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Table 4. The observed data for HD211833. 



Projected Position Hour Normalized Uniform Disk 
MJD Baseline Angle Angle V 2 Ang. Size 

(m) (deg) a (hr) b (mas) 



53173. 


400 


277. 


.40 


114.8 


-2. 


.85 





.1831 ±0. 


.0422 


1 


.219 ±0. 


074 


53173. 


.424 


286. 


.57 


106.1 


-2. 


26 


0. 


1034 ±0 


.0223 


1 


327 ±0. 


053 


53173 


.455 


296. 


33 


96.0 


-1 


.53 


0. 


.1006 ±0 


.0219 


1 


289 ±0. 


051 


53173. 


478 


302. 


.55 


88.5 


-0 


.97 





.1254 ±0. 


.0259 


1 


.214 ±0. 


051 


53175 


400 


279. 


81 


112.5 


-2. 


70 


0. 


1858 ±0 


0503 


1 


.205 ±0. 


088 


53175 


443 


294. 


41 


98.2 


-1 


69 


0. 


1640 ±0 


.0405 


1 


180 ±0. 


072 


53175. 


.484 


305. 


.19 


84.7 


-0 


.68 


0. 


.1068 ±0 


.0251 


1 


.239 ±0. 


055 


53176. 


.406 


282. 


.88 


109.6 


-2 


.50 


0. 


.2647 ±0 


.0802 


1 


.077 ±0. 


117 


53176 


.446 


296. 


27 


96.1 


-1 


.54 


0. 


.1610 ±0. 


.0481 


1 


177 ±0. 


087 


53176. 


498 


308. 


.17 


79.6 


-0 


.29 


0. 


.0782 ± 


.0228 


1 


.290 ±0. 


060 


53177. 


470 


303. 


.18 


87.6 


-0 


.90 





.1270 ±0. 


.0546 


1 


.209 ±0. 


117 


53182. 


449 


289. 


23 


35.9 


-1 


07 


0. 


1772 ±0 


0771 


1 


179 ±0. 


143 


53182. 


.485 


297. 


.53 


26.8 


-0. 


.20 





.0790 ± 0. 


0329 


1 


.334 ± 0. 


095 


53183. 


331 


240. 


.18 


64.4 


-3 


.85 


0. 


.2535 ±0. 


.0648 


1 


.286 ±0. 


.112 


53183. 


.375 


263. 


.39 


53.4 


-2. 


.78 


0. 


.1380 ±0 


.0346 


1 


369 ±0. 


075 


53184 


.364 


259. 


.29 


55.5 


-2. 


.99 


0. 


2736 ± 


.0565 


1 


162 ±0. 


086 


53184 


401 


275. 


.39 


46.4 


-2 


.10 





.2227 ±0 


.0466 


1 


.166 ±0. 


074 


53184. 


438 


287. 


68 


37.3 


-1 


21 


0. 


1897 ±0 


0396 


1 


166 ±0. 


065 


53185 


393 


273. 


55 


47.5 


-2 


21 


0. 


1718 ±0 


0548 


1 


.256 ±0. 


104 


53185. 


.425 


284. 


49 


39.9 


-1 


.47 


0. 


.1960 ±0. 


0629 


1 


.169 ±0. 


107 


53185. 


.453 


292. 


24 


33.0 


-0 


.79 


0. 


1874 ±0 


.0590 


1 


151 ±0. 


100 


53185. 


.485 


298. 


95 


24.8 


-0. 


01 


0. 


1621 ±0 


.0512 


1 


165 ±0. 


092 



a PA is east of north. 

b Errors have been normalized as discussed in §4. 
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The uncertainty in the calibrator visibility represents one of the fundamental limitations of 
the system visibility accuracy. However, our current selection of calibrator is sufficiently 
small in diameter that there are no concerns about a varying system calibration due to 
unaccounted-for calibrator surface morphology. 

Second, SED fitting provides us with an accurate characterization of the stellar bolo- 
metric flux. In the case of our check star, HD 211833, the combination of that flux and an 
actual measure of the star's angular size allows for a direct calculation of the star's effective 
temperature in §4 (rather than the model value used to numerically fit the template). Such 
an analysis will also be applied in §5.2 to our primary target, Alderamin, but as will be 
discussed in that section, a single effective temperature is insufficient to characterize the star 
at the level of detail at which we will be examining it. 

We note that the indicated spectral type and luminosity class for HD 211833 is that of 
K3III (Bidelman 1957; Schmitt 1971), or of a K1III (Wright et al. 2003), but that in Table 6, 
we indicate its best fit spectral template was for a K0-1II star with a SED fitting chi-squared 
per degree of freedom of x 2 /DOF= 0.69, indicating perhaps some uncertainty in not just its 
spectra type but true luminosity class. Fits of the 'normal' Kl through K3 giant models from 
Pickles (1998) indicate x 2 /DOF values of 7.01, 6.18, and 9.31, respectfully; the metal-weak 
and metal-poor templates show no obvious improvement over these values. However, as will 
be examined in the next section, luminosity class uncertainty does not impact our analysis, 
in that HD 211833 appears spherical regardless of luminosity class. 



4. HD 211833 - A 'Round' Check Star 

Before we examine our Alderamin data in detail, we will examine the visibility data 
for our check star, HD 211833. HD 211833 is located in close proximity on the sky to the 



Table 5. Results from spectral energy distribution (SED) fits. 





SED 




x 2 


Fbol 


Model Parameters 


Star 


Template a 


xl 


/DOF 


(10~ 8 erg cm" 2 s" 1 ) 


9.SED (mas) b A v 


HD 197373 


F5IV 


1.08 


19.4 / 18 


10.51 ±0.17 


0.412 ±0.019 0.06 ±0.018 


HD 203280 


A7V 


0.90 


61.8 / 69 


258.3 ± 1.42 


1.36 ±0.034 0±0 


HD 211833 


K0-1II 


0.69 


8.3 / 12 


26.04 ±1.67 


1.32 ±0.29 0.44 ±0.042 



a From Pickles (1998). 

b Estimated angular size from SED fitting as described in §3. 
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HD 197373 Net SED Model 



2.0 x 10" 
1.8x10" 
1.6x10" 
^ 1.4x10"' 
- 1.2x10" 
" g 1.0x10"' 
go 8.0 x 10" : 
^ 6.0 x 10" : 
4.0 x 10" : 
2.0 x 10" : 
0.0 x 10° 



1 1 

Input Photometry 
Model Photometry 
Model Net (Source + Trans) SED 



Fit chi-squared PDF: 1.08 (19.4 / 18 DF) 
1.051e-07 +/- 1.71e-09 erg cnT 
Model T eff : 6562 +/- 150 K 
Est Diam : 0.412 +/- 0.019 mas 




1 1.5 
Wavelength (|J.m) 

Fig. 1. — Spectral energy distribution fitting for our calibrator star, HD197373. In the 
upper panel, the vertical bars on the data points are the errors associated with those data 
points; the horizontal bars represent the bandpass of the data point. In the lower panel, 
the fractional residuals (difference between data point & fit point, normalized by that data 
point) are given for each of those data points. 
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calibrator and primary target (A6> < 7°), and the data collected on all three objects was done 
in a contemporaneous fashion (At < l h ), with little change expected in the point-response 
of the instrument from scan to scan. One of HD 21 1833 's primary attributes that led to 
its selection as our check star was a known, low rate of apparent rotation, with v sin % = 2 
km/s (de Medeiros & Mayor 1999). From that low v sini we inferred that the object would 
have a circular appearance upon the sky, due to either an intrinsic low rotation rate, or a 
pole-on viewing aspect. This circular symmetry is independent of uncertainty in the star's 
luminosity class. 

Once normalized values for V 2 have been obtained as described in §2, the simplest 
interpretation is to fit a uniform disk (UD) angular size to the individual V 2 data points 
following Equation 1. For our check star HD 211833, a single UD fit to the 22 V 2 data points 
results in an indicated angular size of 9 UD = 1.250 ±0.009 mas, with a chi-square per degree 
of freedom fitting value of x 2 /DOF=13.90. Examination of the HD 211833 UD data as an 
ellipsoidal sky projection (as will be detailed in §5.1 for Alderamin) results in fit values of 
2a = 1.252 ± 0.030 mas, 2b = 1.172 ± 0.066 mas, and a = 7.8 ± 6.9 deg, but a reduced 
chi-squared of x 2 /DOF=13.68 - which represents both no significant improvement in fit, and 
more importantly, a negligible detection of asymmetry. 

However, given the known rotational velocities of the calibration and check sources in 
this investigation, it is entirely reasonable to expect that examination of our check star as a 
uniform disk as a function of baseline projection angle should result in a x 2 /DOF of 1.0. For 
the sake of this investigation, we will suggest that the true measurement uncertainty of the 
CHARA Array in the utilized operating mode is not fully characterized by merely tracking 
the measurement scatter as discussed in §2, and that the actual error bars should be a factor 
of 3.42 larger than indicated by that scatter. In doing so, the x 2 /DOF for HD 211833's 
uniform disk fit becomes 1.0, and the indicated uniform disk angular size is 1.235 ± 0.015 

Table 6. Stellar parameters for the check star HD211833 as derived from CHARA angular 

size. 



Parameter 



Value 



Units 



Source 



Apparent rotation velocity (usini) 
Spectral Type 



K3III 
K1III 
K0-1II 



2 



km/s 



dc Medeiros & Mayor (1999) 
Bidelman (1957); Schmitt (1971) 



Wright et al. (2003) 



§3 



Parallax (71-) 

Bolometric flux {Fbol) 

Angular size (9) 

Effective temperature (Teff) 

Linear radius (R) 



4.73 ±0.54 
26.04 ± 1.67 
1.235 ±0.015 



4750 ± 80 
28.1 ±3.2 



K 



mas 



erg cm 



mas 



(Perryman et al. 1997) 
This work 
This work 
This work 
This work 
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mas; as before, the ellipsoidal fit does not indicate a statistically significant improvement. 
We will employ this scaling factor for the errors for examination of the Alderamin data in 
light of what should be the appropriate modeling context. The uniform disk angular sizes 
as a function of baseline projection angle are seen in Figure 2. 

From HD 211833's parallax of tt = 4.73 ± 0.54 (Perryman et al. 1997), we may derive 
a linear radius of R = 28.1 ± 3.2i? . The bolometric flux from §3, in conjunction with the 
angular size, may be used to derive an effective temperature of Teff — 4750 ± 80 K (see 
§5.2 for details on this procedure). Both the radius and temperature numbers are consistent 
with our use of a K0-1II SED fitting template in §3. The full characterization of HD 211833 
that results from this ancillary investigation can be found in Table 6. 



5. Alderamin 

5.1. Apparent Stellar Disk 

The normalized values of V 2 for Alderamin for each observation are listed in Table 3, 
with their associated epoch, projected interferometer baseline length, position angle and 
observation hour angle. A V 2 value is given for each observation, representing the individual 
visibility value derived from the two simultaneous measurements made on light output from 
either side of the beam recombination optic. Also given in Table 3 is an angular size for each 
individual V 2 from uniform disk fit (see Equation 1), which for the purposes of this section 
alone will be used to provide an initial evaluation of the data. Some of the V 2 data points 
lie in the non-monotonic region of a uniform disk visibility curve (where V 2 <~ 0.02 leads 
to multiple possible values of 8ud)', we shall assume for this first look at our data that the 
appropriate values of 9 UD lie on the central lobe of a uniform disk visibility function, noting 
that this assumption does not carry over into our physically more appropriate analysis of 
§6. From those values, fitting a single global value of 9 UD to the V 2 data ensemble results in 
a mean uniform disk size of 1.607 ± 0.032 mas with a chi-squared per degree of freedom of 
X 2 /DOF= 4.484. As is readily evident from Figure 3, a position angle- independent fit would 
clearly be poor. 

This discrepancy can be explored by relaxing the assumption of spherical symmetry 
and including the position angle of the observations in the fit. A spherical gaseous star will 
deform when rotating; such a shape projected onto the sky will appear, to first order, as an 
ellipse. For given physical situations, the true geometry of a rotating star will depart from 
that of an ellipsoid at the 5-20% level, and we will return to this in a much more physically 
appropriate way in §6. However, such a fit is useful as a mathematical construct to initially 
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Fig. 2. — Uniform disk (UD) fits for the individual V 2 data points for HD 211833, as a 
function of baseline projection angle. The UD error bars in the figure are derived from V 2 
errors using the scaling described in §4. The square points are an ellipsoidal fit to the data, 
which for HD 211833 is indistinguishable from a straight line. 
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Fig. 3. — As Figure 2, but for Alderamin. For Alderamin, the projection-angle dependent 
ellipsoidal fit is significantly better than a single angular size fit to all of the visibility data, 
as described in 85.1. 
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establish the position angle dependence of our angular size data. Using the basic equation 
for an ellipse, 

6uD{a) = I 2 ■ V , lta2f ^ (2) 

y or sin (a — aoj + o cos^(o; — «oJ 

we may solve for a projection angle-dependent angular size, where 2a and 26 are the major 
and minor axes of the ellipse on the sky in mas, respectively, and a is the orientation angle 
of the ellipse on the sky with a — corresponding to the major axis pointing to the north 
on the sky and increasing to the east of north. Fitting Equation 2 to the data in Table 2, 
we find that 2a = 1.625 ± 0.056 mas, 2b = 1.355 ± 0.080 mas, and a = 3 ± 14deg with 
X 2 /DOF= 1.08 - this is a substantial improvement over the circular fit. An illustration of 
this fit and the UD data is seen in Figure 4. 

Other potential causes for Alderamin's departure from circularly symmetric V 2 data 
may be ruled out. If Alderamin were either a true or line-of-sight close binary star, our 
interpretation of the V 2 variations with baseline length and position angle would be incorrect. 
If a nearby binary were present in the interferometer beam, variations in the instrument's V 2 
would be present in the data set, but as a function of time, and not just baseline projection 
angle. 

We also consider two other potential deviations of the apparent disk of Alderamin from 
that of a uniform brightness distribution. The first, limb darkening, will affect a star's 
observed visibility curve and potentially bias our results. Second, for a rapidly rotating 
star, this phenomenon takes on an additional latitude dependence, often referred to in the 
literature as gravity-darkening (eg. Claret (2000)). As first shown by von Zeipel (1924a), 
the polar zones of stars distorted by rapid rotation will be hotter than their equatorial zones, 
because the poles are closer to the center of the star. The consequential non-uniform flux 
distribution over the stellar surface affects a star's visibility curve. Our expectation is that 
the ellipsoidal fit in this section can be improved upon by accounting for these effects, which 
we will do in §6. 



5.2. Effective Temperature 

Although we may compute a single effective temperature from our data on Alderamin, 
it must be stressed that this will be nothing more than a mathematical construct derived 
from geometrical considerations for the purposes of characterizing the gross properties of 
the star. Rewriting the stellar effective temperature equation in terms of angular diameter 
and bolometric flux Frol, a value of Trff was calculated from the flux and mean Rosseland 
diameter Or, corresponding to the level in the atmosphere where the Rosseland mean opacity 
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Fig. 4. — Data points along the limb of Alderamin for a simple ellipsoidal fit. The dotted 
line is a circular fit for the same data. 
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is unity, using 

where the units of Fbol are 10~ 8 erg/cm 2 s, and 9r, cir, bn are in mas. The error in Teff is 
calculated from the usual propagation of errors applied to Equation 3. The resultant mean 
T E ff for Alderamin is determined here from the flux value given in §3 and the angular 
size data given in §5.1 to be 7700 ± 170K. Previously, a value of 7773K was estimated by 
Gray et al. (2003), which agrees well with our measure, noting again that this value for 
effective temperature is solely derived from geometric considerations and is an inadequate 
true characterization of a stellar surface over which the temperature, in fact, is latitude 
dependent. 

We should note, however, that this determination of Teff differs from the values for 
pole and equator local effective temperatures that result from the Roche fitting in §6. This is 
due to the accuracy with which the overall bolometric flux can be determined for Alderamin, 
using data across the spectrum from the U band (0.3 /im) to longwards of the M-band (5 
/im). In contrast to that, the photometric fitting portion of the approach detailed in §6 
that constrains the pole temperatures of the models is limited by the accuracy with which 
the V and i^-band brightness of Alderamin has been determined, which is 2.44 ± 0.05 and 
1.96 ± 0.05 (Johnson et al. 1966; Cutri et al. 2003), respectfully. 

A larger implication of this result is the potential inadequacy of effective temperatures 
derived from angular diameters at single projections across the disks of rotationally distorted 
stars. As we will see in the next section, this effect can be much more significant than limb 
darkening in ascertaining a star's Teff, an effect which is expected to be routinely considered 
in all studies of stellar effective temperature. 

6. Rapid Rotator Fitting 

The key to understanding the peculiar diameter results for Alderamin lies in its rapid 
rotation. The force of centrifugal acceleration at the equator, resulting from the rotation, 
offsets the effect of gravitation owing to the mass of the star. Under the conditions of 
hydrostatic equilibrium, uniform rotation, and a point mass gravitational potential, we may 
derive the equatorial rotational velocity, assuming we view the star at an inclination angle 
i. As developed in the work by Collins (1963, 1965) and presented in Jordahl (1972), the 



-19- 



equation of shape for such a star under rotation may be written as 

GM GM 1 2jy . Q 2 . 2 



+ -w 2 R(e,wywi 2 e. (4) 



Rp(u) R(6,w) 2 

From Equation 4, we can arrive at an expression for the colatitude-dependent stellar radius 
at a rotation speed u: 



r(0, u) = — - — - - cos 
u sin 



cos 1 (—usin9) + Arc 



(5) 



where u is the dimensionless rotation speed 

8 GM 
27R p 3 (uj) 



^ = u2 ^TTTr^ (6) 



and r(9,u) is the radius normalized to the stellar polar radius for a given u. It is worth 
noting that, in contrast to our elliptical approximation in §5.1, this approach solves for the 
expected shape of the stellar limb using an approach based upon physics rather than merely 
geometry. 

To interpret our interferometric data, we used a Monte Carlo approach which began by 
constructing models of Alderamin based upon rotation u and polar radius R p {oj), sufficient 
to map the entire surface as a function of stellar colatitude and longitude. Model surfaces 
were constructed for the star at intervals of 0.8° in both colatitude and longitude across 
the whole volume. Flux for a given surface area was then computed using the appropriate 
influence of gravity darkening (von Zeipel 1924a; Claret 2003), with T E ff oc g 13 . For the 
models in question, pole temperature T po i e and gravity darkening coefficient (3 were the free 
parameters that characterized this effect, following the relationship between local effective 
temperature, f3, and local effective surface gravity: 

Ti OC ai(0) = T po i e f - — J (7) 

\9poleJ 

as detailed in Collins (1965) and Domiciano de Souza et al. (2002). These models were then 
mapped onto the sky, through the use of two additional free parameters describing orienta- 
tion, inclination i and on-sky rotational orientation a, with limb-darkening appropriate for 
these model stars as indicated by Claret & Hauschildt (2003) applied at this point as well. 

For comparison with interferometric data, the 2-dimensional model star projected onto 
the sky was Fourier transformed to provide model V 2 's for comparison with all observed vis- 
ibility data points simultaneously. This transformation took into account a mild bandwidth 
smearing effect due to the data being taken through a broad K s filter, by repeating the 
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calculation at a number of points through the filter and averaging the results. Additionally, 
the temperature of each area element of a given model's sky projection would be used to 
compute a contribution to the overall apparent flux density from the star in both the V 
and K s bands. The sum total of the apparent flux density was then compared to measured 
V and K s band photometry for the star, thereby providing a constraint upon T po i e for the 
models. 

Thus, for a given set of six randomized free parameters {u, R p (u),i, a, (3, T po i e }, a 101,000 
point volume surface was generated, projected upon the sky, rotated and the resultant image 
Fourier transformed for comparison to each of the observed V 2 data points, and a x 2 /DOF 
calculated. The multidimensional downhill simplex method optimization code from Numer- 
ical Recipes (Press et al. 1992) was then utilized to derive the best {u, R p (u), i, a, (3, T po i e } 
solution from the random starting point, a process that took typically 500 iterations. In 
contrast to the earlier generation reduction code used in van Belle et al. (2001), this analysis 
compares the model and observed data in Fourier space, rather than image space, which will 
result in a more accurate result. In particular, some of the assumptions regarding uniform 
disk geometry found in van Belle et al. (2001) are no longer invoked; the consideration of 
gravity darkening discussed above is possible only with this approach. 

An exhaustive search of the rotating star parameter space was used to explore the 
X 2 /DOF space through optimizations of over 1,000 random starting points. Furthermore, a 
static grid of {u, i} values was explored for optimal {R p (ou), a, (3, T po i e } values to ensure that 
no local minima were trapping the optimization code. The grid consisted of 1,000 points 
spread uniformly over the space enclosed by u — [0 : 1], i — [0 : 90] and was run multiple 
times with random {R p (ui), a, (3, T po i e } seed values, to ensure full mapping of the resultant 
{u, i} x 2 /DOF surface. 

Once our best {u, R p (uj), i, a, (3, T po i e } solution and its associated x 2 /DOF was estab- 
lished, errors for the individual parameters were derived. This was done by exploring 
the confidence region boundary through a modified version of our optimization code that 
searched the parameter space about our optimum six-parameter fit for appropriate increases 
in x 2 /DOF. Each run of this modified code would target one of the six parameters for max- 
imum deviation from its best-fit value, adjusting the six parameters towards that goal while 
maintaining the A% 2 /DOF constraint. As with the original code, the modified code would 
start with randomized seed values of the six parameters, which in this case were slight de- 
viations off of the best fit. Once done, the multidimensional downhill simplex method code 
would iterate to meet the Ax 2 /DOF condition and maximize the target parameter devia- 
tion. Through approximately a thousand runs of the code for each of the six parameters in 
question, the full confidence region boundary was explored and the appropriate error value 



- 21 - 



was established for those parameters. 

Trial runs of the x 2 /DOF minimization technique using artificial data sets from synthetic 
stars were able to fully recover the initial four parameter characterization for the original 
synthetic star. The model data sets covered a wide range of position angles, from 5deg 
to 175 deg in 5deg steps, and with visibility errors slightly better than in the Alderamin 
dataset, which on average are 4% per measurement. 

The x 2 /DOF surface resulting from the Alderamin dataset is plotted in Figure 6, where 
{Rp(u>), a, [3, Tpoi e } are optimized for minimum x 2 /DOF for a given pair of {u, i} coordinates. 
The six parameter best-fit model's appearance upon the sky is plotted in Figure 7. The best- 
fit value of % 2 /DOF=2.16 is slightly higher than the value for the ellipsoidal fit found in §5.1, 
but this is due primarily to the additional complication of the fit found in incorporating the 
V and K s band photometry constraints. From our solution values for these dimensionless 
parameters, we used values of 2.00 ± O.15M for the mass of Alderamin (derived in §7), and 
a parallax of n = 66.84 ±0.49 mas (Perryman et al. 1997), to extract 'real world' values such 
as rotation velocity and rotational period. 

The difference between the primary and secondary axes of our best fit model is r e (uj) — 
r p {uj) = 201 ± 32 microarcseconds, with an oblateness ratio of 1.298 ± 0.051. Our derived 
value for v sini of 283 ± 19 km/s is in reasonable agreement with the larger spectroscopically 
determined values, as presented in §1. Since our technique of mapping the surface (and 
in particular, the limb) of the star is sensitive to the highest velocity material of the star, 
we were not surprised that our v sin i is on the high end of the spectroscopically determined 
numbers; apparent velocities from spectra have to account for rotationally broadened spectral 
lines convolved across the entire surface of the star, and could potentially underestimate 
values because of this approach. The dominant source of error in our technique is the mass 
estimate, which we will discuss in the next section. The linear sizes are well constrained by 
the Hipparcos parallax, which has only a 0.4% quoted error. 

7. Discussion 

The breadth and depth of the parameters presented in Table 8 allow for a detailed 
examination of the present and past state of Alderamin. As discussed in Reiners (2003) and 
references therein, for spectral types later than A2, a gravity darkening value of (3 = 0.08 
is expected for these stars due to their outer convective envelopes. Our best-fit value of 
(3 = 0.078±o!o59 is consistent with this expectation and lends a degree of confidence to the 
best-fit model. Although the errors on our model value are generous, they do indicate that 
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Fig. 5. — Squared visibility as a function of projected baseline angle for our observed data 
and the best fit model. 
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Fig. 6. — x 2 I DOF surface for Alderamin as a function of rotation u and inclination i. 
Contour lines are in steps of A% 2 / DOF = 0.5 up from the minimum of \ 2 / DOF=2.16 
found at {u,i} = {0.9585,88.2}. 
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Fig. 7. — Best fit 3D model of Alderamin projected onto the sky. The polar regions have a 
temperature of ~ 8440^0 K, with equatorial regions being approximately 7600 K. 
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a gravity darkening value for a radiative envelope of f3 — 0.25 is not consistent with the 
best-fit model. (Running the data reduction code found in §6 and restricting (5 to a value of 
0.25, we find a global minimum with x 2 /DOF=5.01.) 

The radii associated with Alderamin are somewhat larger than the typical values ex- 
pected for a main-sequence A7V star. Our previous investigation of the similarly rapidly 
rotating star Altair (van Belle et al. 2001), with an identical spectral type of A7IV-V, in- 
dicates polar and equatorial sizes that are 30-50% smaller. These increased sizes relative 
to Altair indicate to us that Alderamin is slightly evolved, perhaps being more adequately 
classified as an A7IV; this finding is consistent with an age of 0.82 Gyr for Alderamin as 
quoted by Rieke et al. (2005). This same study found only marginal evidence for excess 
flux at 24 jum, and indicates that there is negligible excess at 2.2 jum, and is consistent with 
the 25 /iin null result of Laureijs et al. (2002). These results are evidence that the observa- 
tions presented here examine the photosphere of the star alone and are not contaminated by 
contributions from a circumstellar disk. 

We can compare our data on Alderamin to the models of Girardi et al. (2000), which 
follow stars of a given mass through their evolution, predicting gravity and temperature; from 
the values of log(g) predicted for these models, we may derive linear radii. We can compare 
the location of Alderamin in radius-temperature space to these predicted tracks, as seen 
in Figure 8; solar metallicity tracks were used, given Alderamin's near-solar metallicity of 
[Fe/H]=0.09 (Gray et al. 2003), and the tracks all start an age of log(T)=7.8 and are stepped 
in increments of A log(T)=0.05. From a simple examination of the location of Alderamin on 
this plot, three new aspects of the star appear to be revealed: First, it's evolutionary status 
of tracking off of the main sequence is confirmed; second, the mass of the object appears to 
be ~ 2.00 ± O.15M , which is consistent with Malagnini & Morossi (1990); and third, its 
age appears to be log(T) ~ 8.9, roughly 800 Myr, which is consistent with the finding of 
Rieke et al. (2005). However, basing these interpretations upon these models in particular 
may be suspect, since the impact of Alderamin's extreme rotation probably alters its specific 
isochrone. We may illustrate this by considering the rotation history of the object through 
conservation of angular momentum. 

The moment of inertia for a star may be written as 

/ = k 2 MR 2 (8) 

where k 2 is the radius of gyration (St§pieh 2000); k 2 = 0.20 for the fully convective case, 
whereas k 2 = 0.05 for the fully radiative case. Although Alderamin's value for k 2 is not 
known in detail, we may consider it for the moment to be constant over the star's recent 
evolution off of the main sequence. From the Girardi et al. (2000) plots in Figure 8, we 
may estimate the average linear size of Alderamin to have increased from ~ 1.6 — 1.8R Q to 
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its current value of ~ 2.54i? . However, by conservation of angular momentum, the star's 
rotation speed when it was this smaller size would have been v e /v c ~ 0.92 — 0.98, very 
nearly rotational breakup speed. Such a previous speed is not impossible from dynamical 
considerations alone, but is far greater than 90% of breakup, which has been argued to be 
the expected upper limit due to star formation considerations - although this limit is not 
borne out by the observations (Stassun et al. 1999; Rebull 2001). 

As such, one of two circumstances may have affected the rotation history of Alderamin, 
independently or in unison. First, the moment of inertia may have changed through changes 
in the radius of gyration as the ratio of convective to radiative portions of the star changed. 
Second, the Girardi et al. (2000) evolutionary tracks potentially do not properly describe 
the radius history of a rapid rotator such as Alderamin. Considerable work on the impact 
of rapid rotation upon the evolution of massive stars (M > 9M ) has been done by Maeder 
& Meynet (eg. Maeder & Meynet (2000)), but ~ 2M solar metallicity model tracks do 
not appear readily available. Alternatively, rotation speeds in excess of 90% of v cr u may be 
allowed in extreme cases such as these. 



8. Conclusions 

We have measured the visibility varying due to the apparent oblateness of Alderamin's 
disk upon the sky and modeled those data with an appropriate Roche model. This approach 
allows for an interferometric measurement of the true stellar rotation velocity and latitude- 
dependent temperature structure, which in turn enabled a more detailed investigation into 
the underlying rotation environment of this star, and its angular momentum history, than 
could be afforded by previous spectroscopic measurements of vsini. The determination of 
the star's gravity darkening from spatial data alone is a unique challenge to stellar models 
and consistent with those theoretical expectations. Such rotational speed and shape de- 
terminations can potentially be also useful for evaluating stellar seismology data (Gizon & 
Solanki 2004). Furthermore, we have demonstrated a technique that can recover a level of 
detail on rapidly rotating stars that heretofore had been out of reach of direct observational 
techniques. Verifying a larger test interferometric cohort around less rapidly spinning stars 
undoubtedly could be key to advances in stellar science. 

As detailed in van Belle et al. (2004), a simple examination of the rotational velocity 
catalog collated by Bernacca & Perinotto (1973) indicates there are over 70 known bright 
(V <4) main sequence stars in the northern hemisphere that are rapid rotators with vsini > 
200 km s -1 ; examination of bright (V <8) evolved objects in de Medeiros & Mayor (1999) 
that have vsini > 15 km s" 1 indicates there are over 70 potential targets as well. Objects 
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Fig. 8. — Alderamin as located on the radius-temperature evolutionary tracks of Girardi 
et al. (2000). Solar metallicity tracks were used, consistent with Alderamin's metallicity of 
[Fe/H]=0.09 (Gray et al. 2003), and the tracks all start an age of log(T)=7.8, being stepped 
in increments of A log(T)=0.05. 



-28- 



that fit these criteria should exhibit apparent flattening of their disks at the rs 10% level. 
Clearly there are numerous opportunities to implement this technique with the new genera- 
tion of long-baseline optical and infrared interferometers such as the CHARA Array, NPOI, 
and VLTI, which all have multiple baselines allowing the required stellar disk projection 
measurements to be made in much shorter observing times. Our CHARA Array follow-up 
observing campaign of other rapidly rotating stars already has initial results that support 
this promising line of research. 

We would like to thank Mel Dyck for first suggesting to us the possibility of utilizing 
oblateness measurements to derive rotational velocity, Doug Gies for a large array of useful 
comments on our manuscript, and Antoine Merand for particularly useful suggestions re- 
garding bandwidth smearing as it pertained to the analysis in §6. Portions of this work were 
performed at the California Institute of Technology under contract with the National Aero- 
nautics and Space Administration. This research has been supported by National Science 
Foundation grants AST-0205297 and AST-0307562. Additional support has been received 
from the Research Program Enhancement program administered by the Vice President for 
Research at Georgia State University. 

9. Appendix - a priori Oblateness Estimation for Rapid Rotators 

Outside of the rigorous mathematical analysis of observed V 2 data for rapid rotators, it is 
useful to have a shorthand approximation of what the expected oblateness for a rapid rotator 
should be. This is particularly useful to developing target lists of these types of objects for 
the CHARA Array and other interferometers. The force of centripetal acceleration at the 
equator, resulting from the rotation, offsets the effect of gravitation owing to the mass of the 
star. Under the conditions of hydrostatic equilibrium, uniform rotation, and a point mass 
gravitational potential, we may derive the equatorial rotational velocity, assuming we view 
the star at an inclination angle i. Under these conditions, we have that 



where Rb and R a are the apparent minor and major stellar radii, and M is the stellar mass (cf. 
Brouwer & Clemence (1961), Elliot & Nicholson (1984), Baron et al. (1989)). Using catalog 
measurements of vsmi (eg. Glebocki & Stawikowski (2000)), and a reasonable estimate for 
M and R (for use as Rb) derived from spectral type, an estimate of the size ratio Rb/R a may 
be established. 




(9) 
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As an example, Malagnini & Morossi (1990) estimate the mass of Alderamin at 1.90 ± 
O.29M ; an A7IV-V star should have an approximate radius of R = 2.1R & (Cox 2000); 
in conjunction with a spectroscopic wsini estimate of ~ 245 km/s (Bernacca & Perinotto 
1970), we find that Rb/R a should be approximately ~ 1.21, which is good agreement with 
the solutions presented in Table 8. 
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Table 7. Alderamin's parameters derived from the data and assembled from the literature. 



Parameter 


Value 


Units 


Reference 


Values from the literature 








Spectral type 


A7IV-V 




(Bidelman 1957) 


Parallax (jr) 


66.84 ±0.49 


mas 


(Perryman et al. 1997) 


Bolometric flux (Fbol) 


2583 ± 14 


10~ 8 erg cm -2 s _1 




Mass (M) 


2.00 ±0.15 


mas 


This work 


Metallicity ([Fe/H]) 


0.09 




(Gray et al. 2003) 


Ellipsoidal fit 








2a R 


1.625 ±0.050 


mas 




2b R 


1.355 ±0.099 


mas 




Position angle (a) 


3± 10 


deg 




Ra 


2.62 ±0.08 


Rq 




Rb 


2.18 ±0.16 


R& 
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Table 8. Alderamin's parameters derived from the gravity- and limb-darkened Roche 

modeling. 



Parameter 



Value 



Units 



Primary Fitting Parameters 



Apparent polar radius (r p (cu)) 


p,7^+°- on9 


mas 


Position angle (a) 


17 2~^'^ 


deg 


Gravity darkening (/?) 


0.0841^ 




Inclination (i) 


oo 9+1.8 
oo.Z ^3 3 


deg 


Polar temperature (T po i e ) 


84401$ 


K 


Dimensionless velocity (u) 


0.9585^ 6 7 





Chi-squared per degree of freedom (x 2 /DOF) 
Derived Values 

Apparent equatorial radius (r e (cj)) 
Polar radius (R p (uj)) 
Equatorial radius (R e (uj)) 
Oblateness (r e (u) /r p {uj)) 
Radii difference (R e (u) — R p (uj)) 
Fractional breakup velocity (v e /v c ) 
Equatorial velocity (v e ) 
Critical velocity (v c ) 
Apparent velocity {v sin 2) 
Period (P) 
Mass (M) 



2.16 

0.876718;^ 
2.175 ±0.046 
2.823 ±0.097 
1.298 ±0.051 
0.649 ±0.107 

n O9O7+0.0482 
-0.0232 

283 ± 19 
342 ± 13 
283 ± 19 
12.11 ±0.26 
2.00 ±0.15 



mas 

Rq 
Rq 

Rq 

km / s 
km / s 
km / s 
hours 

M 



